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Abstract 

This paper makes use of a one-dimensional kinetic model to investigate the nonlinear longitudinal 
dynamics of a long coasting beam propagating through a perfectly conducting circular pipe with 
radius The average axial electric field is expressed as {Ef) = —{d/dz){4>) = —ebgodXb/dz — 
Sb 92 r‘^d^\b/dz^, where go and g 2 are constant geometric factors, Xb{z,t) = / dpzFb{z,pz,t) is 
the line density of beam particles, and Fb{z,pz,t) satisfies the ID Vlasov equation. Detailed 
nonlinear properties of traveling-wave and traveling-pulse (solitons) solutions with time-stationary 
waveform are examined for a wide range of system parameters extending from moderate-amplitudes 
to large-amplitude modulations of the beam charge density. Two classes of solutions for the beam 
distribution function are considered, corresponding to: (a) the nonlinear waterbag distribution, 
where Fb = const, in a bounded region of p^-space; and (b) nonlinear Bernstein-Green-Kruskal 
(BGK)-like solutions, allowing for both trapped and nntrapped particle distributions to interact 
with the self-generated electric field {Ez). 
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I. INTRODUCTION 


High-energy accelerators and transport systems have a wide variety of applica¬ 

tions ranging from basic research in high energy and nuclear physics, to applications such 
as spallation neutron sources, medical physics, and heavy ion fusion. As a consequence, 
it is increasingly important to develop an improved understanding of collective processes 
and the nonlinear dynamics of intense charged particle beam systems. While there has 
been considerable progress in three-dimensional numerical and analytical investigations of 
the nonlinear Vlasov-Maxwell equations describing intense beam propagation, there is also 
considerable interest in the development and application of simplified one-dimensional ki¬ 


netic models to describe the longitudinal nonlinear dynamics of long coasting beams 




15[ | in linear (linac) or large-major-radius ring geometries. The present paper employs the 


one-dimensional kinetic formalism recently developed by Davidson and Startsev [1^ for 
a long coasting beam propagating through a perfecting conducting circular pipe with ra¬ 
dius r^. In Ref, jl^ the average longitudinal electric field is expressed as (Ez) {z,t) = 
— (d/dz) (0) {z,t) = —ebgodXb/dz — ebg 2 rl,d^/dz^, where is the particle charge, go and gz 
are constant geometric factors that depend on the location of the conducting wall and the 
shape of the transverse density prohle, and Xb{z,t) = J dpz Fb (z,pz,t) is the line density. 
In a previous application of the ID kinetic formalism developed in Ref. j^, the analyses in 
Ref. [l5[ assumed that the longitudinal distribution Fb (z,pz,t) corresponded to a so-called 


waterbag distribution 


16l-ll9||. where Fb = const, within moving boundaries in the phase 


space {z,Pz). The weakly nonlinear analysis in Ref. [1^ showed that disturbances moving 


near the sound speed evolve according to the Korteweg-deVries (KdV) equation j20l-l24| . The 


classical KdV equation, which arises in several areas of nonlinear physics in which there are 
cubic dispersive corrections to sound-wave-like signal propagation, also has the appealing 
feature that it’s exactly solvable using inverse scattering techniques. 

While the analysis in Ref. reveals many interesting properties of the nonlinear evo¬ 
lution of longitudinal disturbances in intense charged particle beams, it is restricted to the 
weakly nonlinear regime. In the present analysis, we remove the restriction to the weakly 
nonlinear regime, and make use of the ID kinetic model developed in Ref. iJ], allowing for 
moderate to large-amplitude modulation in the charge density of the beam particles. The 
organization of this paper is the following. In Sec.lTTl the ID kinetic model 1^ is briefly 
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reviewed fSec. lII A|) . and exact (local and nonlocal) nonlinear conservation constraints are 
derived fSec. lIIB]) for the conservation of particle number, momentum, and energy per unit 
length of the beam, making use of the nonlinear Vlasov equation for Fh {z,Pz, t) in Eq. (ITI) . 
and the expression for {E) {z,t) in Eq. ([5]). Removing the assumption of weak nonlinear¬ 
ity made in Ref. 1^, Sec. lIIIl focuses on use of the fully nonlinear kinetic waterbag model 
fSec. lIII A|) to investigate detailed properties of nonlinear pulse-like (soliton) or periodic 
traveling-wave disturbance propagating with constant normalized velocity M = const, rela¬ 
tive to the beam frame fSec. lIUBj) . In normalized variables, Z' = Z — MT and T' = T, the 
waveform of the disturbance is assumed to be time-stationary [d/dT' = 0) in the frame mov¬ 
ing with velocity M = const, relative to the beam frame. Nonlinear solutions are examined 
over a wide range of system parameters, including regimes where the modulation in beam 
line density Xb exceeds 50%, corresponding to a strongly bunched beam. Finally, in Sec.HVl 
we examine the kinetic model based on Eqs. (jH]) and flTU]) [equivalent to Eqs. ([I]) and ([2])] 
for an even broader class of distribution functions Fb {z,pz,t), recognizing that Eqs. ([9]) and 
(ITOl) are Galilean invariant. [Keep in mind that the variables {z,pz, t) are in the beam frame, 
where the particle motion is assumed to be nonrelativistic.] Introducing the appropriately 
scaled variables (see Sec. lIV|) Z' = Z — MT, 1% = Id, — M, T' = T, where M = const., 
we transform Eqs. ([2]) and fllOl) to primed variables, and look for solutions that are time 
stationary [d/dT' = 0) in the frame moving with velocity M = const, relative to the beam 
frame. The analysis in Sec. lIVI parallels the original Bernstein-Greene-Kruskal (BGK) for¬ 
mulation of BGK solutions to the ID Vlasov-Poisson equations 25|, [2^, except for the fact 
that Eq. ffTIl . which connects the effective potential {</>) [z, t) to the line density Xb [z, t), has 
a very different structure than the ID Poisson’s equation used in the original BGK analysis. 
Depending on the choices of trapped-particle and untrapped-particle distribution functions, 
the kinetic model described in Sec.lIVI supports a broad range of nonlinear pulse-like (soli¬ 
ton) solutions and periodic tranveling-wave solution that have stationary waveform in the 
frame moving with velocity M = const, relative to the beam frame. Similar to Sec. lIII Bl 
the modulation on beam line density can have large amplitude, corresponding to a strong 
bunching of the beam particles. Specihc examples are presented in Sec. lIVI corresponding to 
nonlinear periodic traveling wave solutions. 
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II. THEORETICAL MODEL AND ASSUMPTIONS 


This section provides a brief summary of the one-dimensional kinetic g-factor model 


fSec. lII A]) developed by Davidson and Startsev |l^ to describe the nonlinear longitudinal 
dynamics of a long coasting beam propagating in the z-direction through a circular, perfectly 
conducting pipe with radius r^. The ID kinetic Vlasov equation for the distribution function 
Fb{z,pz,t) is used fSec. lII B]) to derive several important conservation laws (both local and 
global) corresponding to conservation of particle number, momentum, and energy per unit 
length of the charge bunch. The results in Secs. lII Al and HT^ form the basis for the nonlinear 
traveling-wave and traveling-pulse solutions studied in Secs, [nil and HVl 


A. Theoretical Model and Assumptions 


This paper makes use of a one-dimensional kinetic model IJ] that describes the nonlinear 
dynamics of the longitudinal distribution function Fb{z,Pz,t), the average self-generated 
axial electric held (Ez) {z, t), and the line density {z, t) = f dpz Ff, {z,pz, t) , for an intense 
charged particle beam propagating in the z-direction through a circular, perfectly conducting 
pipe with radius r^. For simplicity, the analysis is carried out in the beam frame, where 
the longitudinal particle motion in {z,Pz) phase space is assumed to be nonrelativistic, and 
the beam intensity is assumed to be sufficiently low that the beam edge radius and 
rms radius Rb = -1- have a negligibly small dependence on line density 

Aft. Furthermore, properties such as the number density Uft (r, z, f) of beam particles are 
assumed to be azimuthally symmetric about the beam axis {d/dd = 0), where x = rcosO 
and y = rsin6* are cylindrical polar coordinates. Finally, the axial spatial variation in the 
line density Aft {z, t) = 27r fp dr rUb (r, ;s, t) is assumed to be sufficiently slow that -C 1, 
where d/dz kzLp is the inverse length scale of the z-variation. 

Making use of these assumptions, it can be shown that the one-dimensional kinetic equa¬ 
tion describing the nonlinear evolution of the longitudinal distribution function Fb {z,pz,t) 
and average longitudinal electric held (Ez) {z,t) can be expressed in the beam frame correct 
to order as |l^ 


9 F 
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and 


gfe / 17 \ _ Ubo d , d^Xb , . 

m/ Xkodz Aw ^ ^ 

Here, 6^ and mj, are the charge and rest mass of a beam particle, and Aw = const, is a measnre 
of the characteristic line density of beam particles, e.g., its average valne. Moreover, the 
constants U^q and U ^2 have dimensions of speed-sqnare, and are defined by 


7'/2 _ 


^bogocl 


‘^b2 — 


\0g2cl 


mb ^b 

where go and g 2 are the geometric factors defined by [1^ 


(3) 



In obtaining Eqs. (PP)-®, a perfectly conducting cylindrical wall with (r = r^, z,t) = 0 
has been assumed. 

For purposes of illustration, we consider the class of axisymmetric density profiles 
Hb (r, z, t) of the form 


Ah 


rib = 


Trr; 


^/(^), 0^r<rw 
0, n <r 


( 6 ) 


Here, Xb = f dpz Fb {z,pz,t) = 2 tt dr rub {r, z,t) is the line density, rb is the beam edge 
radius, assumed independent of A;,, and / {r/vb) is the prohle shape function with normaliza¬ 
tion Jq dxxf (x) = 1/2. As an example, for / (r/v h) = {n -|- 1) (1 — r'^/rl^ , n = 0,1, 2, • • •, 
over the interval 0 ^ r < r^, it can be shown that 14 1 


= In ( ^ ) + E 


n+l 


n -|- 1 


^^m{m + n + 1) ’ 


92-^ 


n+l 


l + ^1- lb 

, rb2 J ^^m{m + n + 2)r^ 


(7) 


( 8 ) 


(n + 2) + 

From Fqs. (jS])-®, we note that n = 0 corresponds to a step-function density prohle; 
n = 1 corresponds to a parabolic density prohle; n ^ 2 corresponds to an even more 
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sharply peaked density profile; and that the precise values of go and g 2 exhibit a sensi¬ 
tive dependence on profile shape [^. Finally, for the choice of shape function f {r/r^) = 
{n -|- 1) (1 — r'^/r'lY ', n = 0,1, 2, • • •, it is readily shown that the mean-square beam radius 
is Rl = A^^27r /J’’" drrr'^ni, = {n + 2)~^ r^. 


B. Conservation Relations 


Equations ([T]) and ([2]) possess several important conservation laws, both local and global, 
corresponding to conservation of particle number, momentum, and energy per unit length. 
For present purposes we express (Ez) {z,t) = — (d/dz) {(f)) {z,t). Equations ([I]) and ([2]) then 
describe the evolution of Fb {z,pz,t) and (0) {z,t) according to 


where Vz = Pz/ mb and 


5 P ^ F _ 

at az 


d {(f>) dFb 
dz dpz 


0 , 


( 9 ) 


f) f) 

ei^ (4) = mtU^—Ni, + 


( 10 ) 


Here, 


Nb {z, t) = f dpz Fb {z,pz, t) (11) 

■) 

is a dimensionles measure of the line density Xb {z,t), and A^o = const, is the characteristic 
(e.g. average) value of line density. 

It is convenient to introduce the macroscopic moments 


NbVb = Nb {vz) = Xbo j dpz VzFb , ( 12 ) 

Nb{v:) = Xbo^ I dpzV^Fb, (13) 

where Nb = Xb/Xbo is defined in Eq. (ITT]1 . and Vb{z,t) = (f dpz VzFb) (f dpz Fb) is the average 
axial flow velocity in the beam frame. Note that the effective ’pressure’ Pb {z, t) and ’heat 
flow’ Qb {z, t) are defined (relative to the average flow velocity 14) by 


Pb (z, t) = XboNbrUb ({vz - Vbf) ^ 
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and 


Qb {z, t) = XboNbrrib {{v^ - Vbf) = 

where Vb {z,t) is the average flow velocity deflned in Eq. (fl^ . 

We now make nse of Eqs. ([2]) and flTUD to derive the local and global conservation laws 
corresponding to the conservation of particle nnmber, momentnm, and energy per nnit 
length of the beam. The snbsqnent analysis applies to the two classes of beam systems: 

(a) a very long, finite-length charge bnnch {Lb S> r^) with W {z —>■ ±oo,t) = 0; and (b) a 
circnlating beam in a large-aspect-ratio {Rq 3> r^) ring with periodic bonndary condition 
Nb {z + 27 iRq, t) = Nb {z, t) as the beam circnlates aronnd the ring with major radnis Rq. 
(Here, z can be viewed as the arc length aronnd the perimiter of the ring with large radins 
Ro.) 

Number Conservation: From Eqs. ([2]), IfTT]) and flT^ . operating on Eq. ([2]) with / dp^ • • •, 
and integrating by parts with respect to we obtain 

+ = ( 16 ) 

where Nb{z,t) = \b{z,t) /Am is the normalized line density, and Vb{z,t) is the axial flow 
velocity [Eqs. flTT]) and flT^ ]. Eqnation flTBD is a statement of local nnmber conservation, 
i.e., the time rate of change of the local density, dNb/dt, is equal to minus the derivative 
of the local flux of particles, — {d/dz) {NbVb). If we integrate Eq. flTBl) over applying the 
boundary conditions described earlier in this section, we obtain 

-JdzNb = 0, ( 17 ) 

which corresponds to the global conservation of the number of beam particles. 

Momentum Conservation: We now operate on Eq. ([9]) with A^^ / dpzPz ■ ■ ■, where pz = 
mbVz, and make use of Eqs. ffT^ and ffTH . This gives 

d d d 

■^NbmbVb + -^^brrib {vzVz) + ^bNb-^ (</) = 0 , ( 18 ) 

where — {d/dz) ((</>)) is deflned in Eq. ffTUI) . and we have integrated by parts with respect to 
Pz to obtain Eq. flTSD from Eq. ([9]). Equation flTSD can be expressed in an alternate form by 
making use of Eq. flTUD to eliminate Cb {d/dz) {</>) and combine Eqs. ffT^ - flTT)) to express 
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Ntm.1, {v^v,} = JVi.miViV'i, + A„'Pi,, 


( 19 ) 


where Vb {z, t) is the average flow velocity, and Pb {z, t) is the effective pressure of the beam 
particles. Substituting Eqs. (ITII]) and flT^ into Eq. flTS]) . we obtain 


f) f) ( f)JV (9^ At 

-NbrribVb + — {NbmbVbVb + Xbo^Pb} + NbUib \ + U^rl ' 


dz 


^ b2 w 


dz"^ 


= ^^NbrribVb + ^ [NbmbVbVb + X^Pb + ImbU^Nl 


d 


P^mbUl^rl 


N, 


d^Nb 

dz'^ 


1 

2 


'ONb 

dz 


2 

= 0 . 


( 20 ) 


Note that Eq. 0201) expresses the local force balance equation in the form of a local conser¬ 
vation relation for the momentum density of a beam fluid element. Moreover, integrating 
Eq. fl20l) over 2 ; and applying the boundary conditions described earlier in Sec.HTl gives 


J dz NbmbVb = 0 , ( 21 ) 

which corresponds to global momentum conservation. 

Energy Conservation: We now operate on Eq. (E]) with / dp^ \mbv‘l ■ ■ ■ and make 
use of Eqs. (ITTD - (IT5D and p^ = mbVz- Integrating by parts with respect to p^, we readily 
obtain 


^\Nbmb (vf) + ^\Nbmb {vf) + eb^^NbV = 0 . (22) 

From Eqs. flTII]) and flT^ . some straightforward algebraic manipulation gives 


9{(p) d 

o ^b^b o 

oz oz 

d 




^ <1 mbU'^o {NbNbVb) + 


\,,,d^Nb , dNbdNb 

^^bVb r. n + 


dz"^ 


dz dt 


(23) 


Substituting Eq. fl2^ into Eq. fl22D and rearranging terms, we obtain 

















A 

dz 


A 

di 


+ 


(vl) + ^rribU^QN^ - ^rribU^^rl 


I ^Nbrub (vl) + rubUboN^Vb + mbUb2r 


2 _2 
w 


dz 


\,,,d^Nb , dNbdNb 

^VfcKfc ^ ^ + 


dz"^ 


dz dt 


= 0, (24) 


which corresponds to local conservation of energy. Global energy conservation follows upon 
integrating Eq. fl24p over z, which gives 


^ & jiAfimi. («.?) + = 0. (26) 

Note that Eq. fl25p describes the balance in energy exchange between particle kinetic en¬ 
ergy and electrostatic held energy. Moreover, the hnal two terms in Eq. fl2^ correspond 
to electrostatic held energy, and the term proportional to U^q is positive, whereas the term 
proportional to Ub 2 is manifestly negative. Because of the negative sign of the third term 
in Eq. fl25l) . note that any increase in [dNb/dz^ averaged over z must compensated by a 
corresponding increase in the hrst two terms in Eq. (125 1) . 

To summarize, the local conservation laws in Eqs. (ITHll . (lT8l) and fl241) . and the global 
conservation laws in Eqs. (1171) . fl2T]l and fl25ll . provide powerful nonlinear constrains on the 
evolution of the normalized line density Nb, momentum density NbmbVb, and kinetic en¬ 
ergy density Nbmb{v‘l) /2 = NbmbVb / 2 -|- Xb^Pb/^. Furthermore, these conservation con¬ 
straints are exact consequences of the ID nonlinear Vlasov equation ([2]) for Fb{z,Pz,t), 
where eb{d/dz) (0) {z,t) is dehned in Eq. fITUP . and f/^Q and f /^2 expressed in terms of 
the geometric factors go and g 2 in Eq. ([3]). 

Finally, the energy balance equation ([22]), the momentum balance equation fllSp . and the 
continuity equation fITBl) can be combined to give a dynamical equation for the evolution of 
the ehective pressure Pb {z, t) of the beam particles. We make use of Wm?, = NbmbVb + 
Xb^Pbi and express 


Nbmb (vfj = Nbmb ((u^ - Vb + H)^) 

= NbmbVb^ + SNbmbVb ((u, - Vf) + Nbmb ((u, - Vbf) 

= NbmbVb^ + SNbVbXboPb + Qb, (26) 
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where Qb is the effective heat ffow defined in Eq. Without presenting algebraic details, 
some straightforward manipulation of Eq. that make use of Eqs. fll6l) and flTS]) then gives 

| + + + = ( 27 ) 

To summarize, Eqs. (ITBD . (I2U]) and flTH) describe the self-consistent nonlinear of W {z,t), 
Vb{z,t) and Pb{z,t). In the special case where the heat flow contribution {d/dz)Qb is 
negligibly small in Eq. the pressure Pb {z,t) evolves approximately according to 


Pb 


3n?!^ = o, 


I d d 

[m ^ ^ dz 

The continuity equation (ITB]) can be expressed as 


( 28 ) 




Nb + Nb 


m 

dz 


= 0 . 


Combining Eqs. fl27|l and 


we obtain 


( 29 ) 


'd 

dt^^^dz. 




= 0 , 


( 30 ) 


which can be integrated to give the triple-adiabatic pressure relation [Pb/N^) = const. 
Therefore, for negligibly small heat flow in Eq. 0271 ) . the macroscopic fluid model obtained 
by taking moments of the ID Vlasov equation ([9]) closes, and the nonlinear evolution of W, 
Vb and Pb is described by Eqs. OTB]) . (I2UD and (1501) . 

In Sec. mil we discuss a particular choice of distribution function Fb {z,pz,t), correspond¬ 
ing to the so-called waterbag distribution, for which the heat flow Qb {z, t) is exactly zero 
during the nonlinear evolution of the system. In this case, the closure is exact, and the 
nonlinear evolution of the system is fully described by Eqs. (ITO]) . (1201) and (1501) . 


III. COHERENT SONLINEAR STRUCTURES OBTAINED FROM KINETIC 
WATERBAG MODEL 

The ID kinetic g-factor model based on Eqs. ([I]) and ([2]) can be used to determine the 
nonlinear evolution of the beam distribution function Fb {z,pzi t) for a broad range of system 
parameters and initial distribution functions. In this section, we examine Eqs. ([T]) and 
([2]) for the class of exact solutions for Fb{z,Pz,t) corresponding to the so-called waterbag 
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distribution in which Ff, {z, pz, t) has uniform density in phase space fSec. lIII Al) . The subclass 
of coherent nonlinear traveling-wave and traveling-pulse solutions with undistorted waveform 
are then examined fSec. lIUBj) for disturbances traveling in the longitudinal direction with 
constant normalized velocity M = const. 


A. Kinetic Waterbag Model 


Equations ([I]) and ([2]), or equivalently, Eqs. ([2]) and flTUD constitute the starting point in 
the present ID kinetic description of the longitudinal nonlinear dynamics of a long coasting 
beam. The detailed wave excitations associated with Eqs. ([2]) and flTU]) of course depend 
on the form of the distribution function Ff,{z,Pz,t). For small-amplitude perturbations, 
Eqs. ([1]) and ([2]) support solutions corresponding to sound-wave-like disturbances with signal 
speed depending on Ubo and the momentum spread of Fb, and cubic dispersive modihcations 


depending on Ub 2 


li. 


In this section, we specialize to the class of exact nonlinear solutions for Fb {z,Pz,t) to 
Eq. ([I]) corresponding to the waterbag distribution |l5l4l9| 


1 A = const. —nibV^ (z, t) < pz < nibV^ {z, t) , 

(31) 

0, otherwise , 

for —cx) < < cxo (long coasting beam in linear geometry) or 0 < < 2'kRq (large-aspect- 

ratio ring with major radius Rq). In Eq. 0311) . the distribution function Fb = A remains 
constant within the boundary curves —and +mbVi^, and zero outside. The boundary 
curves —mbV,R {z,t) and +mbV,^ {z,t), assumed single-valued, of course distort nonlinearly 
as the system evolves according to Eqs. ([I]) and ([2]) [or equivalently, Eqs. ([2]) and flTU]) ]. We 
integrate across the two boundary curves in Eq.[3l]by operating on Eq. ([T]) with 

rmbV~{l+e) 

lim / dpzPz---, cind lim / dpzPz---, (32) 

e—>-0^" (1—e) Jm^V^{1—e) 

where pz = rubVz. Integrating by parts with respect to and taking the limit e ^ O’*", we 
obtain for the nonlinear evolution of the boundary curves (z, t) and (z, t) 




( 33 ) 
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(34) 


—K+ -—(E) 


where {Ez) is defined in Eq. (|2]). 

For the choice of waterbag distribution in Eq. fl3ip . we calculate several macroscopic fluid 
quantities [see also Eqs. flTTD - flT^ ] corresponding to line density 


Afe = y dpz Fb = Arub ) , (35) 

axial flow velocity 

H = J dpz VzFb = \ (e+ + Vb-) , (36) 

beam particle pressure 

Pb = mbj dpz {vz - Vbf Fb = -^rubA ( 14 + - 14 -)^ = 
and beam particle heat flow 


Qb = mb J dpz {vz - Vbf Fb = Q. (38) 

Note that the heat flow is exactly Qb = 0 for the choice of waterbag distribution in Eq. ([21]). 

Making use of the dynamical equations for and {z,t) in Eqs. (122]) and (1221) . 

where {Ez) is defined in Eq. 02]), some straightforward algebra shows that \b{z,t), Vb{z,t) 
and Pb {z, t) evolve according to 




d.. 1 dPt , (Ul,d , , Ulrla^\ 




m dz i Afeo dz Am dz^ I 


Note from Eqs. (I27|) and (I2T]) that Pb (z, t) can be expressed as 


Pb{z,t) = ^A^ {z,t) , 


(39) 

(40) 

(41) 


(42) 
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where Pbo = const, and Xbo = const, represent the characteristic (e.g., average) valne of the 
pressnre and line density of the beam particles respectively, and Pbo/Xl^ = 1/12 (mbA)'^ = 
const ., where A is the constant phase-space density in Eq. 0311) . By virtne of the fact that the 
heat flow Qb {z, t) = 0 exactly for the choice of distribntion fnnction Fb {z, p^, t) in Eq. it 
is not snrprising that Eqs. are identical to the macroscopic flnid eqnations flTBD . 020]) 

and (I5U]) . obtained in Sec. IllB1 where Eq. (I5UD has made the assumption of negligible heat 
flow in Eq. fl27)) . Note here that Nb {z, t) and Xb {z, t) are related by Nb {z, t) = Xb {z, t) /Am- 
For present purposes, we introduce the effective thermal speed UbT associated with the 
waterbag distribution in Eq. fl3Tl) dehned by 


Jj2 _ 3Ao 

Xb,mb ’ 

and the nomalized (dimensionless) fluid quantities p {z, t) and U {z, t) dehned by 


(43) 


p = Nb -I 


Xb — X 


bO 


X 


u = 


Vb 


bO 


+ U^t) 


1/2 


(44) 


1 /o 

In Eq. fl44]l . (f/^Q -|- U'^rp) ' is the effective sound speed associated with the geometric factor 
go and the thermal speed UbT- Furthermore, we introduce the scaled (dimensionless) time 
variable T and spatial variable Z dehned by 


Making use of the macroscopic equations fl39l) . (HOI) and (02]), and the dehnitions in Eqs. (031) - 
(0^ . it is straightforward to show that the continuity equation (IHHj) and force balance equa¬ 
tion (001) reduce in dimensionless variables exactly to 




(46) 


U? 


bT 


d d (1 2 1 

^ ^ ( 2 ^ + + 2 Uio + U^T 


2 ^ 

p + ITTpjP 


= 0 . 


(47) 


dZ'‘ 

The huid description in scaled variables provided by Eqs. (0f)1) and (071) is exactly equivalent 
to the kinetic description provided by Eqs. ([T]) and ([2|) for the choice of waterbag distribution 
in Eq. (0T]) . 
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B. Coherent Nonlinear Traveling Ware and Pulse Solutions 


Within the context of the present ID model, Eqs. and fH7|) can be used to investigate 
detailed properties of collective excitations over a wide range of system parameters. For 
example, in the weakly nonlinear regime, for small-amplitude disturbances moving near the 
sound speed (f/^g -|- ^ ; Eqs. and fHTD can be shown to reduce to the Korteweg- 

deVries equation |^, which exhibits the generation and interaction of coherent structures 
(solitons) for a wide range of initial density pertuibations rj {Z,T = 0) 7 ^ 0 2^. While the 
analysis in Ref Q has several interesting features, the results are limited to the weakly 
nonlinear regime where <C 1 and \U\ -C 1. 

In this paper, we examine Eqs. 0461) and fl47|) in circumstances where there are not a priori 
restrictions to small amplitude, i.e ., rj = {Xh — Am) /Am is allowed to be of order unity, as 
long as Xb/Xho > 0, which corresponds to r] > — 1 . Furthermore, we look for solutions to 
Eqs. 046 p and 047p that depend on Z and T exclusively through the variables Z' = Z — MT 
and T' = T, where M = const, is the normalized pulse speed measured in units of the 
sound speed (t/^g -|- Making use of d/dZ = djdZ' and djdT = d/dT' — Md/dZ' 

and looking for time-stationary solutions {d/dT' = 0) in the frame of reference moving with 
normalized velocity M = const. , Eqs. 046 h and 0471) for r] {Z') and U {Z') become 


d 


[(-M + U)p + U] = t), 


( 48 ) 


d 


2 


MU + r] + - 


U? 


bT 


Integrating with respect to Z', Eqs. 


2 Ui, + 

and 04^ give 


2 

V H- - 

I ^ QZ'^ 


= 0 , 


( 49 ) 


— Mr] + {1 + r])U = const ., 


( 50 ) 


-U'‘ - MU + + 


'-'bT 


2 


-rf + 


bT 


d‘^rj 


= const. 


(51) 


which relate rj (Z') and U (Z'), where Z' = Z — MT . 

The solutions for p (Z') and U (Z') to Eqs. flSOp and 0^ depend on the values of the 
constants in Eqs. fISUp and 0441) . For present purposes we consider boundary conditions such 
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that U = Q when r; = 0, and r]'' = 0 when U = Q and r; = 0. In this case the values of the 
constants in Eqs. flHJ]) and fIFI]) are zero, which gives 


U = M- 


T] 


1 +r] 


(62) 


d‘^r] 


+ 1-{U - Mf - -M^ + 7] + - , 


v"} = o 


Substituting Eq. into Eq. (IS5D . we obtain 


(53) 


l2 1(1 + 1,)" 




r/a 

^bT 


V =0, 


(54) 


2Ul, + Ulj. , 

which is a second-order nonlinear differential equation for the perturbation in line density 
Tj (Z') = [Ab (Z') — Abo] /Abo, where Z' = Z — MT. Some straightforward algebraic manipu¬ 
lation shows that Eq. fIM]) can be expressed in the equivalent form 


d^r] _ d 


(55) 


where V {rj) is the effective potential dehned by 


1 77^ 

VM = -- 


21 + 7 ] 

and the dimensionless parameter e^, dehned by 


{er?7^ -F (1 -1- er) ?? + (l - } , 


(56) 


6 

^ 3+1 + Uij.) ' 


(57) 


is a measure of the longitudinal thermal speed of the beam particle. 

Equations (IH^ and (1561) can be used to determine the solutions for rj (Z') for a broad range 
of dimensionless parameters and M. Furthermore, Eqs. fl5^ and fl5BD have been obtained 
from Eqs. (15Ui) and fIFID for the special class of boundary conditions where 1/ = 0 and 
7]" = 0 when r; = 0 [see discussion prior to Eqs. fl5^ and fl5^ ]. Indeed, we will show below 
that Eqs. (I5H]) and fl5BD support two classses of solutions consistent with these boundary 
conditions. These correspond to: (a) localized (pulse-like) soliton solutions when > 1, 
satisfying rj [Z' = ±cx)) = 0, 17 iZ' = ±cx)) = 0, and [c^^h/^-^^^]z'=±oo ~ (t*) nonlinear 

periodic traveling-wave solutions when < 1, with t] (Z') = rj {Z' + L), and rj {Z' = 0) = 0, 
U (Z' = 0) = 0, and [d‘^7]/dZ''^]^,^^ = 0. 
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In general, the effective potential V (rj) in Eq. flHTD can be expressed as 


where 


^ (^) = [(^ - 77+) (r; - r; ) 


2 1 + 77 


( 58 ) 


77 


± 



± 



(59) 


In Eqs. fl5H]) and (1591) . is restricted to the range 0 < < 1/3, and can satisfy > 1 

or < 1. Examination of Eq. (l59l) shows that 


77 < — 1 , 77 + > — 1 , 


(60) 


for all allowed values of and Furthermore, it’s also clear from Eq. that 


[ 77 + > 0 , for > 1 , 

[^ 77 + < 0 , for < 1 . 


(61) 


Recall that rj = {Xb — Am) /Am- Then A^/Am ^ 0 implies that 77 ^ —1 is the region of interest 
physically for solutions to Eq. ([T]). Note that Eq. (|55ll has the form of a dynamical equation 
of motion, with 77 playing the role of displacement, Z' playing the role of time, and V ( 77 ) 
playing the role of effective potential. Multiplying Eq. 0551) by drj/dZ' and integrating, we 
obtain 


\ (^) + ^ (^) = ^ = ( 62 ) 

Equation 062p plays the role of an energy conservation constraint, and can be integrated 
to determine 77 (Z') for the effective potential V {rf) defined in Eq. 058 p . We now examine 
solutions to Eq. 062p for the two cases indentified earlier: < 1 and —1 < 77 + < 0; and 

> 1 and 77 + > 0. 

Nonlinear Traveling-wave Solutions < 1 and — 1 < 77 + < 0).' Figure [T] shows a 
schematic plot of V ( 77 ) versus 77 for the case where < 1 and —1 < 77 + < 0. For 
purpose of illustration, the values of the specific parameters in Fig. [T] have been chosen to 
be = 0.09 and ct = 4/15 in plotting V ( 77 ) versus rj. The corresponding values of 77 +, 
77m and V (rjm) are 77 + = —0.882, rjm = —0.715, and V ( 77 ^) = 0.126. For different choices 
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of values for er and < 1, the shape of the V (rj) versus rj curve is qualitatively similar 
to that shown in Fig.[TJ Referring to Fig.dl when the effective energy E (the red horizonal 
line in Fig.[T] lies in the interval Q < E < V {r]m), Eq. 0661) supports nonlinear periodic 
solutions for rj [Z') that oscillate as a function of Z'. Here V {rjm,) is the local maximum of 
V (p), which occurs dX rj = rjm in Fig.[TJ Depending on system parameters, these nonlinear 
traveling-wave solutions can have large amplitude, representing a siginihcant modulation in 
beam line density. 

Referring to the discussion preceding Eq. flSSl) . the boundary conditions used to derive 
Eqs. fl^ and 0^ from Eqs. 0^ and 0^ correspond to 

ri{0)=0 = ri” (0) (63) 

for the class of nonlinear periodic wave solutions obtained from Eq. when < 1 and 
— 1 < < 0. Furthermore, from Fig. [Hand Eq. flF^ . we note that R (p = 0) = 0 and the 

effective energy E can be expressed as 

E = \w(a)f- (64) 

Typical numerical solutions for r] {Z'), obtained by integrating Eq. with V [r]) specihed 
by Eq. fISS]) . are illustrated in Figs. [2J12] for several values of < 1 and e^, and different 
values of effective energy level E. These correspond to: = 0.36, er = 0, E = 0.005 

and E = 0.0110707 (Fig.[2]); = 0.36, er = 4/15, E = 0.005 and E = 0.020 (Fig.Ej); 

m2 = 0.09, eT = 0, E = 0.005 and E = 0.054483 (Fig.jl); = 0.09, er = 4/15, E = 0.005 
and E = 0.125 (Fig.[5l). Close examination of Figs. [2ll5] shows several interesting trends. 
First, for smaller values of M^, the potential wells are deeper and broader (compare Figs.EK 
and Hti, and Figs.jSt'' and [5^); and for smaller values of e^, the potential wells are deeper 
(compare Figs. [2^ and [3^, and Figs.jl^ and [5^). Furthermore, the nonlinear wave amplitude 
tends to be larger for smaller values of M^ (compare Figs. [2^ and [3^, and Figs.jl^ and [3^), 
whereas the wavelength dependance on M^ and ex tends to be relatively weak (compare 
Figs.lH, El m [5]). In any case, for M^ < 1, it is clear from Figs.[T]l5] that Eqs. (1621) and 
(1331) support a broad class of nonlinear traveling-wave solutions for the theoretical model 
developed here, based on the ID kinetic waterbag model for intense beam propagation. 
Indeed, the modulation of the beam line density is about ±50% for the system parameter 
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in Figs. St and St- 

Nonlinear Traveling-pulse (Soliton) Solutions > 1 and ri~^ > 0).' We now consider 
Eqs. (16^ and Eq. (j58ll [ or equivalently, Eq. (j56|) ] in circumstences where > 1 and ri~^ > 0. 
In this case, the effective potential has the qualitative shape illustrated in Fig. SI which has 
been plotted for the choice of parameters = 9 and er = 1/30. The physically allowed, 
localized pulse solutions (soliton solutions) corresponds to the energy level 

E = 0, (65) 

which is the red horizontal line in Fig. El and boundary conditions 

r] {Z' = ±cx)) = 0 = rj" {Z' = ±cxo) (66) 

discussed prior to Eq. (l58l) . Referring to Fig. El when Eq. (l62l) is integrated forward from 
Z' = —oo where p = 0, the perturbed line density, rj increases monotonically through 
positive values to a maximum amplitude (the soliton amplitude) and then decreases 
monotonically to r/ = 0 when Z' = +cx). The regime where > 1 by a sufficiently large 
amount corresponds to a strongly nonlinear regime where the density compression is large 
with r]^ > 1. On the other hand, when — 1 = e is small with 0 < e -C 1, the soliton 
amplitude is correspondingly small. This will become apparent from the numerical solutions 
to Eqs. and (ESI) consistent with Eqs. and (E 6 ]) presented later in this section in 
Figs.lTlfTUl 

Typical numerical solutions to Eqs. fE^ and (ESI) , subject to Eqs. fE^ and fEB]l . are 
illustrated in Figs.lTlfTOlfor several values of > 1 and er- These correspond to: = 4 

and er = 0 (Fig.[7]); = 4 and ct = 4/15 (Fig.E]); = 1.2 and ct = 0 (Fig.ED; and 

= 1.2 and = 4/15 iFig. fTUl) . Close examination of Figs.lTHTUl shows that the soliton 
amplitude increases with increasing (compare Figs. [7] and E] with Figs.Eland fTU]) . reaching 
a highly nonlinear regime with r]~^ = 3.0 in Fig.[7]and = 1.735 in Fig.El where = 4. In 
contrast, the soliton width tends show a relatively weak dependence on longitudinal velocity 
spread, as measured by er (compare Fig.[H]with Fig. [71 and Fig. [101 with Fig. [9]). It’s clear 
from Figs.lTlfTni that the soliton solutions to Eqs. fETI) and fEHl) exhibit a strong nonlinear 
dependence on M^, and can correspond to highly compressed line density for sufficiently 
large M^. 
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In the special circumstances where exceeds 1 by a small amount, i.e., = 1 + A 

where 0 < A 1. it is readily shown that Eq. flM|l can be approximated for small rj by 


d‘^ri 


+ 


~M + —€t 


7 ] — I 7] = 0 


(67) 


Equation fl^ can be solved exactly for rj (Z') = Af, (Z') /Abo — 1 to give 


rj (Z') = 




M^ + er) 
Note that the soliton amplitude in Eq. 


sech^ 


-(m^ 

9 V 


1/2 

1) (Z - MT) 


( 68 ) 


is small for = 1 + A with A -C 1. Also, the 


sech^ {■ ■ ■ } pulse shape in Eq. (16 8 h is similar to the soliton pulse shape obtained from the 
Korleweq-deVries equation in the weakly nonlinear regime 15l |. 

Finally, it should be noted that the oscillatory solutions obtained from Eqs. (1581) and 
(1621) when > 1 and the energy level E in Fig. El is negative with Vmin < E < 0 are 
not considered here. These solutions are unphysical because they oscillate about a positive 
non-zero average value of i] = Afe/Afeo — 1 > 0, rather than oscillate about f] = Ab/Abo — l ~ 0, 
as occurs in Figs. [2}E] when < 1. 


IV. COHERENT NONLINEAR STRUCTURES OBTAINED FROM FULLY KI¬ 
NETIC G-FACTOR MODEL 

The kinetic waterbag model developed in Sec. lIIIl of this paper has clearly demonstrated 
the rich variety of coherent nonlinear structures supported by the ID kinetic model based 
on Eqs. ([9]) and flTOl) [or equivalently, Eqs. ([I]) and ([2])] for the specihc choice of waterbag 
distribution Eb {z,pz, t) in Eq. (16T1) . In this Section, we examine solutions to Eqs. (E]) and (ITO]) 
for an even broader class of distribution functions Eb {z,Pz,t), recognizing that Eqs. ([6]) and 
fllOp are Galilean invariant. That is, if we transform variables to a frames of reference moving 
with constant longitudinal velocity Vq = const, according to z' = z — Vot, p'^ = Pz — fnbVo, 
t' = t, then in the new dynamical variables {z',p'^,t'), the equations for Eb {z',p'^,t') and 
(0) {z',t') are identical in form to Eqs. ([9]) and (ITOj) . Time-stationary solutions {d/df = 0) 
in the new variables {z',p'^,t') then correspond to undistorted traveling-wave or traveling- 
pulse solutions moving with constant velocity Vq = const, in the original variables {z,Pz,t). 
The present analysis of Eqs. ([6]) and flTUD parallels the original Bernstein-Greene-Kruskal 


19 









(BGK) formulation of BGK solutions to the ID Vlasov-Poisson equations 2^, 12^, except 


for the fact that Eq. ffHJ]) . which connects (0) {z,t) to the line density Xb has a very 

different structure than the ID Poisson equation. 

Referring to Eqs. ([2]) and flTU]) . we introduce the scaled dimensionless variables (Z, P^jT) 
dehned by 


V ( uio + uiT \~^ z ^ ( uio + uiA Ub2t 

p = _ ^ _=_ — _= y ('69') 

+ (Ul + USrP " 

where U^q and U ^2 dehned in Eq. ([3]), and U^rp = const, is the longitudinal velocity 
spread characteristic of the distribution function Fb. We futher introduce the dimensionless 
distribution function Fb {Z, P^, T) dehned by 


Fb 





(70) 


where Am = const, is the characteristic line density of the beam particles, e.g., the average 
value. From Eqs. fl69ll . flTHl) and the dehnition of line density Xb = J dpzFb, it follows that 
the perturbation in line density rj = Xb/Xbo — 1 can be expressed as 


T] {Z, T)= I dPz Fb {Z, Pz, T) - 1, (71) 

where the Pz integration covers the range — oo < P^ < cx) in Eq. fl71l) . Transforming variables 
according to Eqs. fl39]) and flTU]) . and making use of Eq. flTT]) . it is readily shown that Eqs. ([3]) 
and flTOll can be expressed in the new variables as 


and 


dT ^ ^ dZ dZ dPz 


(72) 


Ip = T] + 


d^T] 


where ip {Z, T) is the normalized (dimensionless) potential dehned by 


Ip = 




{Ul + Ulr) 


1/2 


(73) 


(74) 
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Equations dza) and f|73p . wlioro tj and J dPzFb arc rnlatnd by Eq. f|7ip . constituts cou¬ 
pled nonlinear equations describing the self-consistent evolution of the distribution function 
Fh {Z, Pz, T), normalized potential V’ {Z, T), and normalized perturbed line density 7] {Z, T). 
Equations (j7Tl) - (j73jl are fully equivalent to the original dynamical equations (B-dli. and 
can be used to investigate ID kinetic properties of the nonlinear beam dynamics over a wide 
range of system parameters. 

Keeping in mind that Eqs. fl72p and 0731) are Galilean invariant, if we transform Eqs. 072p 
and fITSp from the variables {Z, Pz, T) to a frame moving with normalized velocity M = const. 
according to Z' = Z — MT, = Vz — M, T' = T, then Eqs. fl72p and (1751) have exactly the 
same form in the new variables, with {Z,Pz,T) replaced by {Z', P'^,T'), i.e.. 


and 


^ 

dT' ^ ^ dZ' dZ' dP^ 


(75) 


Ip = 7] + 


d’^rj 


Here, P^ = VP, and Fb {Z', Pp, T') and r] (Z', Pp, T') are related by 


(76) 


r, (Z', T) = J dPl Ft (Z', P'„ U) - 1. (77) 

Therefore, the traveling-pulse or traveling-wave solutions that have stationary prohle shape 
in the primed variables (Z', Pp, T') are determined by setting d/dT' = 0 in Eqs. (j75i) - (j77l) . 
Setting dPb/dT' = 0 in Eq. (1751) gives for Fb (Z', PP) 


yfdPb _ dxp dPb 
^ dZ' dZ' dPp 


(78) 


where (Z') and r] (Z') solve Eq. (1751) . and rj (Z') is related to Fb (Z', PP) by Eq. (ITTp . We 
introduce the energy variable W dehned by 


»" = pvf + (Z') . 

Then the solution to Eq. (1751) for Fb (Z', PP) can be expressed as 


(79) 


Ft (Z\ U) = P (»") e (U) + Ff (W) e (-V ’), 


(80) 
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where 


Note from Eq. (17^ that 


0 (K) 


forV'>0, 

[o, /orf/'< 0 . 


(81) 


dV^ = ±dW'/ [2 iW' - , (82) 

where + corresponds to > 0, and — corresponds to fd' < 0. Snbstitnting Eqs. fISUl) and 
fIST]) into Eq. fl77D gives 


h = 


dW'^ 


F> {W') + F< (hh') 


1 , 


(83) 


[2 {W'- 

which relate the pertnrbation in beam line density rj [Z') to the potential ip {Z') and the 
distribntion fnnctions Pp' {W') and Pp^ (W'). 

Fignre [TTl shows an illnstrative plot of the potential ip {Z') as a fnnction of Z'. Depending 
on the valnes of the energy W and range of Z', there are three classes of particle orbits: (a) 
particles that are reflected from the potential; (b) particles that are trapped and nndergo 
periodic motion; and (c) passing (nntrapped) particles that don’t change direction, bnt pass 
over the potential maximnm, Erst slowing down and then speeding np dnring the motion. 
For the trapped particles and the reflected particles, it follows that P^ (hhO = P^ (^0 
that 


pTr {W) = P<, {W) + {W') = 2p<^ (W') = 2pp {W) (84) 

and 

Fr., (W) = Flf (W) + Flf (W) = 2FS,, (.W) = 2F^; (W) . (86) 

On the other hand, for the passing (untrapped) particles, P^^ [W') and Ppf^ (hhO can be 
specihed independently, depending on whether the particles have VP > 0 or 1/' < 0, respec¬ 
tively. 

The form of ip {Z') shown in Fig.fTTI corresponds to a stationary isolated pulse in primed 
variables, with ip {Z' = ±cx)) = 0. By contrast, Fig.fT^ shows a plot of ip {Z') versus Z' for 
the case where ip {Z') has a periodic nonlinear wave structure with 
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«i(Z' + L)=^(Z'). 


( 86 ) 


From Fig. [121 trapped particles with energy W in the range 


i^niin < Ih"' < V’max (87) 

exhibit periodic motion. On the other hand, passing particles with energy W in the range 
(see Fig.[T2D 


> i^niax ( 88 ) 

correspond to untrapped particles that pass over the potential V’ {Z'), periodically speeding 
up and slowing down, but not changing their direction of motion. Furthermore, for the 
nonlinear periodic waveform for the potential 'ip {Z' + L) = p) (Z') shown in Fig.[T21 it follows 
from Eqs. flTbp and flHdl) that the waveform for the perturbation in line charge also satishes 
rj (Z' + L) = rj (Z'). Here, 7] (Z') is related to ^p (Z') and the trapped-particle and untrapped- 
particle distribution functions by Eq. (15^ . which gives 


1 + r/ 



pTr jW) 

[2 


+ 



Fun (W') 

[2 {W - ■ 


(89) 


In Eq. (l89il . the integration over the trapped-particle distributuion Ftt iW) is over the 
interval of W corresponding to V’mm < < W < ipmax, and the integration over the 

untrapped particle distribution Fun iW) is over the interval of W corresponding to ipmax < 
Ip < W < CX). 

Equations (17^ and flS^ can be used to determine detailed properties of self-consistent 
nonlinear periodic solutions for rj [Z') and ip {Z') for a broad range of choices of Ftx (hFO 
and Fun (hFO- Furthermore, depending on system parameters, the amplitudes of the wave 
perturbations can range from small to moderately large amplitude. For purposes of illustra¬ 
tion the procedure for solving Eqs. (1761) and fl8^ for the case of nonlinear periodic solutions 
for rj (Z') and ip (Z'), we consider the special case where Ftt iW) = 0, and the untrapped 
distribution function has the monoenergetic form 


Fun (W) = (IE' - W[j) , 


(90) 
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where = const ., A = const., and Wlj > 'ipmax (see Fig.[T2]). Substituting iW) = 0 
and Eq. fl^ into Eq. flS^ readily gives 


1 + ,,(Z') = 


A 


(91) 


For present purpose, we choose the normalization constant A in Eq. (l91|) such that the 
line density perturbation ?] (Z') = A/, (Z') /Am — 1 and potential perturbation -i/ (Z') are 
simultaneously zero for all Z', i.e., rj (Z') = 0 for all Z' when -0 (Z') = 0. From Eq. 0911) . 
this readily gives A = 1 for the value of the constant A. Squaring Eq. flHT]) and solving for 
Ip {Z') when A = 1 readily gives 


rP = Wp 


1 - 


(92) 


(1 + >))t 

Note that Eq. (|9^ determines ip [Z') as a function of rj (Z'), which can be substituted into 
Eq. d7H]) to solve for rj (Z'). 

Similar to the analysis in Sec. lIIIBl for the class of nonlinear periodic traveling-wave 
solutions with rj {Z' + L) = rj (Z') and ip (Z' + L) = ip (Z'), we examine Eqs. fITBl) and fl^ for 


the case where the boundary conditions correspond to rj iZ' = 0) = 0 and = 0. 

Substituting Eq. (|^ into Eq. (1761) we readily obtain 


d^ri 


+ rj 


which can also be expressed as 


wp 




d^rj dV 
dZ'"^ dr] 

where 


(93) 


(94) 


drj 


Wlj 


ii + vY\ 


{i + vY 


'rj^ + (2-Wlj)rj+{l 



(95) 


Note that Eq. fIMll has the form of a dynamical equation of motion, with rj playing the role of 
displacement, Z' playing the role of time, and V {rj) playing the role of an effective potential. 
Making use of Eq. fl^^ . it is readily shown that 
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(96) 


ffiV_ 2W'i 

'W~ (Tm)? ’ 

and 




r] + 


l + T] 


- 1 


1 


T] 


{r,+ |l-2<]} , 


(97) 


2 ( 1 +,)= 

where the constant of integration in Eq. fl97p has heen chosen so that E (77 = 0) = 0. 

Close examination of Eqs. shows that Eq. fl^lD snpports oscillatory solntions for 

rj (Z') abont rj = 0 provided [d'^V/> 0, or eqnivalently, 


2iy^ < 1. (98) 

When the ineqnality in Eq. fl^S]) is satished, the plot of V (r/) versns r] has the characteristic 
shape illnstrated in Fig.[T3]for the choice of parameter 2hh^ < 1. Here, V [t]) has a minimnm 
at 77 = 0 , and passed throngh zero at 


7/ = 77+ = - [1 - 2 WIj ] , (99) 

where V {rj = r]~^) = 0 [see Eqs. (1971) and fl99]) ]. Similar to the analysis in Sec. lIIIBl Eq. fl91]) 
can be integrated to give the energy conservation relation (1/2) [drj/dZ']^ + V {rj) = E = 
const, [see also Eq. (1621) ]. where E = (1/2) [t]' (0)]^ is the effective energy level. Referring to 
Fig.Hni Eqs. (199]) and (I97|) snpport nonlinear periodic oscillatory solntions for 77 {Z') for E 
in the range 0 < E < Vm, where Vm = V{r] = rjm) is the local maximnm of V ( 77 ) at 77 = rjm- 
For the choice of dimensionless parameter 2 HC = 1/2 in Fig.[T31 it is readilty shown that 
hm = —0.360 and Vm = V {rj = rjm) = 0.014. 

Recall that the primed variables {Z', V/,T') are related to {Z, 14, T) hj Z' = Z — MT, 
14 = 14 — M and T' = T, where M = const, is the dimensionless velocity of the traveling 
wave relative to to the nnprimed frame. Therefore, for a nonlinear wave that is time station¬ 
ary [d/dT' = 0) in the primed variables, it is reasonable to identify IF^ with H 4 = (1/2) 
for a monoenergetic beam. In this case, we make the identihcation 2Wlj < 1, so the condi¬ 
tion for Eqs. (199]) and (1991) to have nonlinear periodic solntions for 77 {Z') [see Eq. (199]) ] can 
be expressed as 
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M^<1. (100) 

Typical numerical solutions for rj (Z'), obtained by integrating Eq. 0941) with V (rj) specified 
in Eq. 097|) . are illustrated in Figs. fTHfTTI for several choices of < 1 and different values 
of effective energy level B. These correspond to: = 0.5, B = 0.005, = —0.360, 

and V(r]„,) = 0.014 (Fig.[IlD; = 0.5, B = 0.054883, = -0.360, and V (rj^) = 0.014 

(Fig. USD; = 0.09, B = 0.05, = —0.764, and V (rjm) = 0.181 (Fig.[Tn]); and = 0.09, 

B = 0.18, rjm = —0.764, and V {rim) = 0.181 (Fig.dTI). Figures ITilfTTl illustrate several 
interesting trends in the nonlinear periodic wave solutions for ri{Z'). [These should be 
compared with the nonlinear periodic wave solutions in Figs.[2}15] obtained in Sec. lIIIl for 
the kinetic waterbag model.] First, for smaller values of M^, the potential wells are deeper 
and broader (compare Figs.fTTh andlTHk with Figs.fTOh and ITTh) . Furthermore, the nonlinear 
wave amplitudes tend to be large for sufficiently large values of energy level B in the potential 
well (compare Figs.fTHI and ITTl with Figs. dH and USD- 


V. CONCLUSIONS 


In this paper, the ID kinetic model developed in Ref. m was used to describe the non¬ 
linear longitudinal dynamics of intense beam propagation, allowing for moderate-to large- 
amplitude modulation in the charge density of the beam particles. Particular emphasis has 
been placed on investigating detailed properties of nonlinear pulse-like (soliton) and periodic 
traveling-wave disturbances propagating with constant normalized velocity M = const, rel¬ 
ative to the beam frame. The ID kinetic formalism was briefly summarized in Sec. Ill A 


and exact (local and nonlocal) nonlinear conservation constraints were derived in Sec. Ill Bl 
for the conserved particle number, momentum, and energy per unit length of the beam, 
making use of the nonlinear Vlasov equation for B), {z,Pz,t) in Eq. ([T]) and the expression 
for (Bz) {z,t) in Eq. ([2]). Removing the assumption of weak nonlinearity made in Ref. |^, 
Sec. lIIIl made use of the fully nonlinear kinetic waterbag model to investigate detailed prop¬ 
erties of traveling nonlinear disturbances propagation with velocity M = const, relative to 
the beam frame. In normalized variables, Z' = Z — MT and T' = T, the waveform of the 
disturbance was assumed to be time-stationary {d/dT' = 0) in the frame moving with veloc¬ 
ity M = const. Nonlinear solutions were examined over a wide range of system parameters 
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for both traveling-pulse (soliton) and nonlinear traveling wave solutions in which the mod¬ 
ulation in beam density was large-amplitude, corresponding to a strongly bunched beam. 
Finally, in Sec.[IV]we examined the kinetic model based on Eqs. (jH]) and flTU]) [equivalent to 
Eqs. ([I]) and ([2])] for an even broader class of distribution functions {z,Pz, t). The analysis 
in Sec. IIVI parallels the original Bernstein-Greene-Kruskal (BGK) formulation of BGK solu¬ 


tions to the ID Vlasov-Poisson equations 23, [2^, except for the fact that Eq. fllOl) . which 


connects the effective potential {(j)) {z,t) to the line density Xb{z,t), has a very different 
structure than the ID Poisson’s equation used in the original BGK analysis. Depending 
on the choices of trapped-particle and untrapped-particle distribution functions, the kinetic 
model described in Sec. HV] supports a broad range of nonlinear pulse-like (soliton) solutions 
and periodic traveling-wave solutions that have stationary waveform in a frame of reference 
moving with velocity M = const, relative to the beam frame. Similar to Sec. lIIIl the modu¬ 
lation in beam line density can have large amplitude, corresponding to a strong bunching of 
the beam particles. Specihc examples were considered in Sec. IIVI corresponding to nonlinear 
periodic traveling wave solutions of Eqs. ([9]) and flTOll . 
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FIG. 1. Illustrative plot of V{r]) verses r] obtained from Eq. (56) for = 0.09 and er = 4/15. 
Here, r]~^ = —0.882, rjm = —0.715 and V{r]m) = 0.126. 
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FIG. 2. For = 0.36, et = 0, r/+ = —0.64, = —0.476 and V{rim) = 0.0355, plots are shown 

for (a) V{i]) verses rj; (b) ri{Z') verses Z' for ri'{0) = 0.1 and E = (l/2)[r/'(0)]^ = 0.005; and (c) 
r]{Z') verses Z' for r/'(0) = 0.2 and E = (l/2)[r/'(0)]^ = 0.02. 
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FIG. 3. For = 0.36, er = 0.8, = —0.442, r/m = —0.306 and V{r]m) = 0.0111, plots are 

shown for (a) V{rii) verses rj; (b) 'r]{Z') verses Z' for V(O) = 0.1 and E = (l/2)[r/'(0)]^ = 0.005; and 
(c) ri{Z') verses Z' for ^(0) = 0.1488 and E = (l/2)[r/'(0)]^ = 0.0110707. 
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FIG. 4. For = 0.09, e-r = 0, = —0.91, = —0.764 and V{r]m) = 0.181, plots are shown for 

(a) V{r]) verses r/; (b) r]{Z') verses Z' for r]'{0) = 0.1 and E = (l/2)[77'(0)]^ = 0.005; and (c) r]{Z') 
verses Z' for r/'(0) = 0.5 and E = (1/2)[//'(O)]^ = 0.125. 


32 









-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 






Z' 


FIG. 5. For = 0.09, er = 0.8, = —0.767, r/m = —0.557 and V{r]m) = 0.0545, plots are 

shown for (a) V{rii) verses rj] (b) 'r]{Z') verses Z' for r]'{0) = 0.1 and E = (l/2)[r7'(0)]^ = 0.005; and 
(c) ri{Z') verses Z' for r]'{Q) = 0.3301 and E = (l/2)[?7'(0)]^ = 0.054483. 
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FIG. 6. Illustrative plot of F(r/) verses rj obtained from Eq. (56) for = 9 and ex = 1/50. Here, 
r]~^ = 6.908, and the energy level E = 0 corresponds to soliton solutions with maximum amplitude 
r/+ = 6.908. 






FIG. 7. Plots of (a) V (r/) verses r/; and (b) r]{Z') verses Z', obtained from Eqs. (56) and (62) for 
= 4, eT = 0 and E = 0, corresponding to soliton amplitudes = 3.0. 
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FIG. 8. Plots of (a) V (r/) verses r/; and (b) r]{Z') verses Z', obtained from Eqs. (56) and (62) for 
m2 = 4, eT = 4/15 and E = 0, corresponding to soliton amplitudes = 1.735. 
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FIG. 9. Plots of (a) V{rf) verses r/; and (b) r](Z') verses Z', obtained from Eqs. (56) and (62) for 
m 2 = 1.2, e-r = 0 and E = 0, corresponding to soliton amplitudes r]~^ = 0.4. 


36 











FIG. 10. Plots of (a) V{rj) verses rj; and (b) r]{Z') verses Z', obtained from Eqs. (56) and (62) for 
m2 = 1.2, ex = 4/15 and E = 0, corresponding to soliton amplitudes r]~^ = 0.297. 
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FIG. 11. Illustrative plot of the effective potential verses Z' occurring in Eq. (33) showing 

the three classes of particle orbits corresponding to (a) passing (untrapped) particles with energy 
IF 3 , (b) reflected particles with energy VF 2 , and (c) reflected or trapped particles (depending on 
the range of Z') with energy W[. The form of 'ijj{Z') in Fig. 11 corresponds to an isolated pulse 
with ijj{Z' —)■ ± 00 ) = 0. 



FIG. 12. Illustrative plot of the effective potential 'ijj{Z') verses Z' for the case where V’(-^0 
a nonlinear periodic waveform with 'ip(Z' + L) = 'il){Z'), where L is the periodicity length. In 
the figure, passing particles with energy W' > ipmax are untrapped, whereas particles with energy 
'ipmin < < tpmax are trapped and exhibit periodic motion in the potential il^{Z'). 
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FIG. 13. Plot of V{r]) verses rj obtained from Eq. (97) for 2Wl^ = 0.5 and = —0.5. Here, 
rjm = —0.36, and V{r]m) = 0.014. Nonlinear periodic solutions for ri{Z') exist for energy level E in 
the range 0 < E < V{r]m)- 
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FIG. 14. Plots are shown for (a) V{r]) verses rj, and (b) verses Z', obtained from Eq. (94) for 

m 2 = 0.5, r/'(0) = 0.1, E = 1/2[7/'(0)]2 = 0.005, t]+ = -0.5, rjm = -0.36 and V{r]m) = 0.014. 
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FIG. 15. Plots are shown for (a) V{rj) verses rj, and (b) r](Z') verses Z', obtained from Eq. (94) for 
m2 = 0.5, ? 7 '( 0 ) = 0.1863, E = l/2[r]'{0)]^ = 0.054883, r/+ = -0.5, rjm = -0.36 and V{r]m) = 0.014. 
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FIG. 16. Plots are shown for (a) V{r]) verses rj, and (b) r](Z') verses Z', obtained from Eq. (94) for 
m 2 = 0.09, r/'(0) = 0.1, E = l/2[r/'(0)]2 = 0.05, r]+ = -0.91, rjm = -0.764 and = 0.181. 
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FIG. 17. Plots are shown for (a) V{rj) verses rj, and (b) rj(Z') verses Z', obtained from Eq. (94) for 
m 2 = 0.09, r/'(0) = 0.6, E = l/2[r/'(0)]2 = 0.18, 77 + = -0.91, rjm = -0.764 and V{r]m) = 0.181. 
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